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Knowledge of atmospheric turbidity coefficients is very important in meteorology, climatology, 
atmospheric pollution monitoring, and in the prediction of solar energy availability under cloudless 
skies. This paper provides a thorough review on a number of atmospheric turbidity indices and on the 
several methods that have been developed in the past few decades. The Lambert-Bouguer-Beer law is 
the basic relationship underlying the derivation of various turbidity indices. Turbidity can be referred to 
as monochromatic (narrow band) wavelengths, broadband wavelengths, and the total spectrum. 
Narrow band turbidity is measured using sun photometers, while the other two turbidities are 
measured using pyrheliometers with broadband pass filters. The Angstrém’s turbidity coefficients 
and the Linke’s turbidity factor are among the most frequently used atmospheric turbidity coefficients. 
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1. Introduction 


Most energy produced in the fusion furnace of the sun is 
radially transmitted as electromagnetic radiation, which is popu- 
larly called sunshine or solar energy [1]. The sun has a tempera- 
ture of 5760 K, and it emits the total power of approximately 
3.8 x 107° W or 6.25 x 10!! Wm”. Since the distance between 
the sun and the earth is large, the amount of solar radiation that 
reaches the outside of the earth’s atmosphere is quite low at only 
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1367 W m~’, which is called the total solar irradiance (TSI) or the 
solar constant [2]. The solar constant is measured on the surface 
perpendicular to the sun’s rays at the average sun-earth distance 
above the atmosphere [3]. The intensity of solar radiation above 
the atmosphere approximately varies sinusoidally over the year 
with amplitude close to 3.3% of the solar constant and it achieves 
its maximum around the first day of January. This variation arises 
from the variation in the distance of the earth from the sun [4]. 
Variation in solar radiation is also caused by the difference in the 
emission intensity from the sun itself [3]. 

Not all solar radiation received at the periphery of the atmo- 
sphere reaches the Earth’s surface [5]. Solar radiation is partially 
depleted and attenuated as it traverses the atmospheric layers [6]. 
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Nomenclature 

B the Angstrém’s turbidity coefficient 

a the Angstrém wavelength exponent 

Tad spectral aerosol optical depth at the wavelength 2 

Vis visibility in kilometers 

Iq diffuse solar radiation (W m~?) 

Ip direct solar radiation (W m~?) 

Ip the irradiance of the extraterrestrial solar radiation 

Isc extraterrestrial solar irradiance at normal incidence 
(1367 W m~?) 

Te threshold irradiance 

Ibn normal direct solar radiation at the earth’s surface 
(W m~?) 

Eo correction term for the earth-sun distance 

bn normal incidence solar radiation at the bottom of a 
dust-free atmosphere 

T5 the aerosol optical thickness at 0.5 um 

Be Angstrom turbidity coefficient obtained using Guey- 
mard method (dimensionless). 

m, relative optical air mass at a given pressure 
(dimensionless) 

Tpin Linke turbidity factor obtained using Pinazo et al. 
method (dimensionless) 

Kab the ratio between the diffuse and direct beam normal 
irradiance 

Zaff the effective wavelength corresponding to each of the 
filters 

Tadd the aerosol optical thickness at a specified 
spectral band 

T; Linke’s turbidity factor (dimensionless) 

Tim monthly average values of Linke’s turbidity factor 

ôr(m,) the integral Rayleigh optical thickness 


Or(2) the relative optical thickness relating to Rayleigh 
scattering by the gaseous molecules 


Op(Z) the relative optical thickness associated with aerosol 
extinction and gaseous absorption other than ozone 

Orx(m,;) the integral Rayleigh optical thickness computed by 
the formula of Kasten [63] 

Tix(m,) Linke’s turbidity factor determined by Kasten’s algo- 
rithm [63] 

Tia(m;) Linke’s turbidity factor determined by Kasten’s algo- 
rithm [68] 

Ora(m,) integral Rayleigh optical thickness computed by the 
formula of [68] 

This Linke’s turbidity factor in the visible region of the 
solar spectrum 

Tan direct hourly normal solar radiation 

Tup Linke turbidity factor obtained using Pinazo et al. 
method (dimensionless) 

TR the mean Rayleigh extinction 

Tw the selective absorption coefficient for atmospheric 
gases 

Tu Unsworth-Monteith turbidity coefficient 

Tram2 Linke’s turbidity factor for an air mass equal to 2 

w precipitable water thickness (cm). 

Kab the standardised value (to zero altitude and total 
column ozone equal to 0.3434 atm cm) of the ratio 
between diffuse and direct irradiance 

B Schiiepp’s turbidity coefficient 

m* an absolute air mass 

CDA the number of clean-dry atmospheres 

y solar elevation angle 

Ti illuminance’s turbidity factor 

Ty Volz’s turbidity factor 

TR Rayleigh scattering coefficient for air moleculares at 
the wavelength 2 

1H absorption coefficient for ozone at 1 

Tk Kastrov’s turbidity coefficient 

T; rational’s turbidity factor 


The presence of aerosol particles in the atmosphere is one of the 
most important factors that affect the amount of solar radiation 
that reaches the Earth’s surface under cloudless skies [7]. Aerosols 
are small particles that are either in the solid or liquid state [8]. 
Common aerosol particles have radius ranging from 1 to 10° nm, 
while very small particles (called Aitken particles) range from 1 to 
10 nm and large particles range from 10° to 10? nm. Particles in 
the range of 10? to 10° nm are called giant particles [9,10]. Based 
on the formation of aerosol dispersions, aerosol particles have 
two categories, namely, dispersed and condensed aerosols [8]. 
Aerosols have two dominant layers in the atmosphere. The first 
dominant layer is near the Earth’s surface at 0-3 km, which is 
affected by natural dust storms and man-made inputs to the 
atmosphere, whereas the second dominant layer is the strato- 
spheric dust layer at 15-25 km above sea level, which is caused 
by volcanic action and cosmic sources [11]. 

Aerosols have direct radiative forcing because they scatter and 
absorb solar and infrared radiations in the atmosphere. Aerosols 
also alter warm, ice, and mixed-phase cloud formation processes 
by increasing droplet concentrations and ice particle concentra- 
tions, consequently causing indirect radiative forcing associated 
with the changes in the microphysical and optical properties of 
clouds. Furthermore, aerosols are highly absorbent of solar radia- 
tion, such as black carbon, and they may reduce cloud cover and 
liquid water contents by heating the cloud and the environment 
where the cloud forms. This is known as the semi-direct effect 
because it is the result of the direct interaction of aerosols with 


radiation that indirectly influences the climate by altering clouds. 
The net effect of aerosols on climate depends on the summation 
of all the three mechanisms above [12]. An atmosphere that 
contains aerosols is described as turbid or hazy. An aerosol-laden 
atmosphere has the property called atmospheric turbidity, which 
depletes direct solar radiation [10]. 

Studies on atmospheric turbidity are important in meteorology, 
climatology, and in atmospheric pollution monitoring. Atmospheric 
turbidity parameters are necessary for determining the amount of 
spectral global irradiance and for designing photovoltaic cells and 
selective absorbers for spectral thermal collectors [13]. Aerosols 
make the atmosphere turbid, and thus, the presence of aerosols is 
commonly quantified by an index of turbidity. In the last century, 
some indices of turbidity were proposed and several techniques 
were developed to determine their values [14]. Among these indices 
are the Angstrém’s and the Schiiepp’s turbidity coefficients, which 
are spectral in nature, and the Unsworth-Monteith’s coefficient 
and Link’s turbidity factor, which are broadband in nature [15]. 
This study mainly aims to investigate atmospheric turbidity 
and then describe available techniques for calculating various 
turbidity indices. 


2. Atmospheric aerosol characteristics 


The word aerosol was introduced more than 80 years ago as an 
analogy to the term hydrosol, a stable liquid suspension of solid 
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particles. Aerosols are products of a complicated totality of 
chemical and physical processes [16,17]. Aerosols are defined in 
their simplest forms as small solid or liquid particles that remain 
suspended in the air and follow the motion of the air within 
certain broad limits. An atmosphere that contains aerosols is 
described as turbid or hazy. An aerosol-laden atmosphere has the 
property called atmospheric turbidity, which depletes direct solar 
radiation [16,10]. The presence of aerosols in the Earth’s atmo- 
sphere has important effects on the transmission of solar radia- 
tion and on the radiative heat transfers in the atmosphere. This is 
due to the fact that aerosols absorb and scatter solar radiation as 
it passes through the atmosphere [14]. Scattering refers to 
processes in which photons change direction after an interaction, 
while absorption occurs when photons are removed from a beam 
light and their energy is converted to an excitation of atoms or 
molecules [18]. Aerosols are physically divided into the following 
three categories according to particle size: 


1. Aitken nucleus with particle radius smaller than 10? nm 

2. Large nucleus with particle radius larger than 10? nm and 
smaller than 10? nm 

3. Giant nucleus with particle radius larger than 10? nm 


Based on the formation of the aerosol dispersions, aerosols fall 
into the following two categories [8]: 


1. Dispersed aerosols: the aerosol materials that are either in the 
solid or liquid state are pulverized into a granular state via 
mechanical pulverization or natural weathering, and are then 
suspended in the air because of the increase in wind force. 

2. Condensed aerosols: the formation of the aerosol particle or 
droplet is caused by the condensation of super-saturated gas 
on the condensation nucleus or by the mixture of different 
gases through photochemical reaction. 


Studies on the physical properties of aerosol particles are 
important for the assessment of the effects aerosol particles on 
the climate and for the development of more accurate remote- 
sensing procedures from satellite sensors [19]. Aerosol properties 
vary with time and from one region to another because of their 
short lifespan and highly temporal and spatial variability mea- 
surements [20,21]. The spatial-temporal variability of aerosol 
characteristics is large and the observational data are so frag- 
mentary that the total content of various types of aerosols is 
difficult to estimate, while the available power estimates of the 
global sources of natural and anthropogenic aerosols are rather 
approximate [17]. 

Kondratyev et al. [17] and Giannakaki et al. [22] reported that 
the types of aerosols depend on their sources and are roughly 
categorized into anthropogenic or natural. Anthropogenic sources 
are determined by human activities, such as industrial wastes 
from chimneys, toxic exhausts from cars, fires, explosions, soil 
erosion in agriculture, and open mining. The estimates of anthro- 
pogenic aerosols are more reliable compared with those of 
natural aerosols, especially in regions along the World Ocean 
and in the continents that are difficult to access. Atmospheric 
aerosols are classified into the following types of natural aerosols 
based on their composition or sources: 


1. Products of sea spray evaporation. 

2. Mineral dust wind-driven to the atmosphere. 

3. Volcanic aerosols (both directly emitted to the atmosphere 
and formed by gas-to-particle conversion). 

4. Particles of biogenic origin (directly emitted to the atmosphere 
and formed by the condensation of volatile organic compounds, 


such as terpenes, and from the chemical reactions between 
these compounds). 

5. Smokes from biota burning on land. 

6. Products of natural gas-to-particle conversion (e.g., sulfates 
from reduced sulfur on the ocean surface with emissions of 
dimethyl sulfide). 


Important types of anthropogenic aerosols also include the 
following: 


1. Direct industrial emissions of particles (soot, smoke, and road 
dust, among others). 
2. Products of gas-to-particle conversion. 


The amount of aerosols in the atmosphere is sometimes 
measured in terms of the number of particles per cubic centi- 
meter. The number of dust particles is usually greater on land 
than over water and during the drier seasons. Moreover, the 
number of dust particles is lower in cold polar or arctic air and is 
higher in tropical air. The number density exponentially decreases 
with altitudes of up to approximately 5 km and remains constant 
at altitudes of 10-15 km. The number density slightly increases in 
a layer centered around 20km because of the thermal and 
dynamic nature of the atmosphere [10]. Tropospheric aerosol 
particles play an important role in the attenuation of solar 
radiation because they strongly affect the Earth’s radiation bud- 
get. The influence of aerosols on radiation that passes through the 
atmosphere cannot be neglected, especially in urban or indus- 
trialized areas, [20,21,23]. 

Aerosols are a minor constituent of the Earth’s atmosphere, 
but they play an important role in the energy balance of the 
earth-atmosphere system. Aerosols can influence the climate in 
two ways [24]. The magnitude of these influences, however, is 
poorly constrained because of limited knowledge on the pro- 
cesses that control the distributions as well as the physical, 
chemical, and optical properties of aerosols [25]. First, aerosol 
particles can reduce solar energy by absorbing or scattering the 
radiation under cloud-free and cloudy conditions, thereby affect- 
ing the energy balance of the Earth [24]. This direct effect is 
related to the optical properties of aerosol particles, such as single 
scattering albedo, size distribution, size variability with relative 
humidity, complex refractive index, and solubility. The levels of 
scattering and absorption depend on the physical and chemical 
characteristics of aerosol particles. Consequently, aerosols modify 
the Earth’s radiation budget and consequently influence the 
warming or cooling of the planet. Second, aerosols alter the 
formation and the precipitation efficiency of liquid-water, ice, 
and mixed-phase clouds, thereby causing indirect radiative for- 
cing associated with these changes in the cloud properties. The 
quantification of indirect radiative forcing caused by aerosols is 
an especially difficult task [26]. 

Radiative forcing of the climate by aerosols may be compar- 
able in magnitude, but opposite in sign, to that of greenhouse 
gases. Greenhouse gases trap radiation and warm the surface of 
the earth, while aerosols absorb and scatter radiation back to 
space, theoretically resulting in a net cooling effect [24]. This 
variability is largely due to the much shorter atmospheric lifetime 
of aerosols compared with that of greenhouse gases [23]. Turbid- 
ity observations had been made since 1923, but their measure- 
ments are inherently difficult and expensive even making some of 
them not possible because of the strong influence of other atmo- 
spheric components [27,15,28]. The accurate determination of 
turbidity normally requires clear-sky spectral radiation data 
obtained using sun photometers or spectroradiometers. These 
instruments are expensive and scarce, and thus, turbidity is 
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instead generally estimated from broadband irradiance measure- 
ments [29]. 


3. Atmospheric turbidity indices 


Aerosols have a relationship with the attenuation of solar 
radiation that reaches the Earth’s surface, and thus, a number of 
methods and atmospheric turbidity indices have been introduced 
and developed in the past decades to evaluate atmospheric 
turbidity and to quantify its influence on the direct irradiance 
on the Earth’s surface [14,30,31]. The Angstrom turbidity coeffi- 
cients « and f and the Schüepp turbidity coefficient B, among 
others, are spectral in nature, such that they reference the 
spectral aerosol optical depth to a common wavelength (i.e., 
1000 nm for f and 500nm for B). On the other hand, the 
Unsworth-Monteith turbidity coefficient Ty and the Link turbidity 
factor T, are broadband in nature, i.e., they are estimated by 
reversing or inverting the models that estimate broadband beam 
irradiance upon measurement. The Unsworth-Monteith turbidity 
coefficient is, by definition, equal to the broadband aerosol optical 
depth, which is the sum of aerosol absorption and scattering 
coefficients integrated over the entire depth of the atmosphere 
[15,30]. The definition of each turbidity index is explained in the 
succeeding sections. 


3.1. Angstrém’s turbidity coefficients (« and ß) 


Under clear sky conditions, atmospheric aerosols generally 
contribute to the largest attenuation of solar radiation in the 
visible range. Scattering and absorption effects cause the attenua- 
tion. Since attenuation effects of the scattering and absorption by 
aerosols are difficult to separate, Angstrém suggested a single 
formula that is generally known as the Ångström turbidity 
formula, which can be expressed as follows [10,21,32]: 


Tal = pre aí) 


where T,,; is the spectral aerosol optical depth at the wavelength, 
B is the Ångström turbidity coefficient, ~ is the wavelength 
exponent, and / is the wavelength in micrometer. 

In Eq. (1), 6, which varies from 0.0 to 0.5 or even higher, is an 
index that represents the amount of aerosols in the atmosphere at 
the vertical direction. The wavelength exponent « is related to the 
size distribution of the aerosol particles. Large values of « indicate 
a relatively high ratio of small to large particles, which vary from 
0 to 4. When the aerosol particles are very small in relation to air 
molecules, « should approach 4, whereas for very large particles, 
a should approach 0. Generally, «x is has a value between 0.5 and 
2.5. The average « for the entire solar spectrum has a value 
between 0.9 and 2.0. An acceptable average value for most natural 
atmospheres is «= 1.3 + 0.5. Table 1 lists the combinations of 
parameters f and «, which may be considered for various degrees 
of atmospheric cleanliness [10,33]. 

The Ångström exponent « is a widely used parameter in 
atmospheric sciences that examines the optical properties of 
aerosol particles. Since the early [34,35] and later [36,37] pub- 
lications of Ångström, where the « parameter was mainly used to 


Table 1 
Parameters for various degrees of atmospheric cleanliness [10]. 


Atmosphere B a Visibility (km) 
Clean 0.00 1.30 340 
Clear 0.10 1.30 28 
Turbid 0.20 1.30 11 
Very turbid 0.40 1.30 <5 


describe the spectral behavior of atmospheric extinction and 
transmission, « is now also applied to a variety of similar but 
slightly different optical properties such as the atmosphere 
scattering or backscattering coefficients. The Ångström exponent 
is very popular because of the simplicity of its equation that 
enables the interpolation or extrapolation of aerosol optical 
properties and because of its connection to particle microphysics 
(related to the mean size of aerosols) [38]. 


3.1.1. Techniques for determining the Ångström turbidity coefficients 
(x and f) 

Angstrém turbidity coefficients « and f can be determined 
using a number of techniques. Traditionally, « and f are derived 
from either spectral direct solar radiation measurements or from 
broadband direct solar radiation and precipitable water measure- 
ments [39,40]. 


3.1.1.1. The Volz method. Volz [41] applies Eq. (1) to sun photometric 
measurements performed in two different narrow bands that are 
chosen to minimize the influence of other atmospheric attenuators, 
such as water vapor. In this case, x can be obtained as follows: 


dinta, _ In(=) (2) 
dnd nz 


where 2, and Az represent the effective wavelengths in the 
corresponding narrow bands [21,26]. 

The chosen wavelengths are usually 380 and 500 nm, where 
weather conditions play a very important role in the distribution 
of pollutants in the air [39,28]. At 380 nm, no molecular absorption 
is observed, while at 500 nm, the ozone has weak absorption [10]. 
Eq. (2) indicates that « is the negative of the slope or the negative of 
the first derivative of t,; versus the wavelength in logarithmic scale. 
After « is derived, f can be calculated in the same wavelength 
interval using the values of « and T,,, at a certain 2 [21]. However, f 
alone can be measured with a single-wavelength Volz instrument by 
assuming that «=1.3. Alternately, f can be measured at A=1 um, 
where the exponent is not included in the calculations. At this 
wavelength, only the water vapor has a weak absorption band. This 
method is very accurate and yields the values of « and f simulta- 
neously [10,39,28]. 


3.1.1.2. Direct method. The direct method enables the calculation 
of the Angstrém turbidity coefficients y and f by using the 
following formula obtained from Eq. (1) [42]: 
d(Int,,] A [dt 
dink — Ty 


(3) 


The parameter f can also be determined from Eq. (1). 


3.1.1.3. Linear fitting. Linear fitting depends on the omission of Ta, 
in the strong absorption bands of water vapor and mixed gases, 
and it determines the Angstrom parameters at remaining 
wavelengths. This fact constitutes the main difference with the 
previous (direct) method. The Angstrom parameters are derived 
by applying the least-squares fit, as written in Eq. (1), to the 
remaining wavelengths, which can be rewritten in the following 
linear form [42]: 


Int,, = —a Ind +Inp (4) 


In this approximation, the slope of the straight line yields a, 
while the intercept provides f. Hence, « is the negative of the 
slope (first derivative) of In Tą, versus In 2. This method has often 
been used for determining aerosol characteristics because of its 
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simplicity [21,43]. Since the Angstrom formula does not exactly 
fit the data points, a second-order polynomial fit for In t,; versus 
In 2 data provides excellent agreement with the measured T,,,, 
resulting in order differences of the uncertainty of the measure- 
ments (~0.01-0.02), while a linear fit yields significant differ- 
ences with the measured T,,;. The second-order polynomial fit is 
applied to the t,, values at all wavelengths and is expressed in 
the following form [26]: 


Inty; = %2 Ind? +01 INA+a9 (5) 


Since « is the negative of the first derivative of In tą, versus In /, 
the following expression can be derived from Eq. (4) to quantify 
the spectral dependence of «: 


& = —2% Ind—a, (6) 


3.1.1.4. Horizontal visibility. Routine measurements of f and « 
with a dual-wavelength sun photometer are conducted at a 
number of locations. When such measurements are unavailable, 
the turbidity parameter f may be determined from the 
measurement of visibility (Vis) in the horizontal direction, if 
such a measurement exists. This visibility is also called the 
meteorological range. For visibilities greater than 5km, the 
turbidity parameter can be determined using the following 
equation developed by McClatchey, Selby, and others [10]: 


B =0.55%(3.912/Vis—0.01162)[0.02472(Vis—5) + 1.132] (7) 


where Vis is in kilometers. In order to predict f}, the value of a 
must first be known. Unfortunately, even with the help of proper 
instruments, the measured visibility is slightly subjective in 
nature. Therefore, the relationship above only yields a gross 
estimate of f. 


3.1.1.5. Rangarajan and Mani method. Rangarajan and Mani [44] 
described a new, simpler, and less expensive method of evaluating 
p from direct and diffused solar radiation measurements. Their 
method has been a routine in principal radiation stations. 


Iq = 0.5Sa + 0.7554 


r = 0S0 SD Pa 

Sa = 1-g(f)™ (9) 
0.606m* 

Se= 6434 me au 

g(B)=1-1.375ß for f < 0.08 (11) 

g(B) = 0.9843-1.1429f for B>0.08 (12) 


where m* is an absolute air mass, J, is the direct solar radiation, 
and l4 is the diffused solar radiation. 


3.1.1.6. Louche et al. method. Louche et al. [39] described a method for 
determining the Angstrém’s turbidity coefficient in a Mediterranean 
location from the total pyrheliometric measurements in the para- 
meterization of the direct solar total irradiance by Bird and Hulstrom 
[45,46], which can be expressed as follows: 


1 C 
P= nD nà) a3) 
where m, is the relative air mass at actual pressure (dimensionless) 


and 


A= Ipn/ (0.975EolscTrToTgTw) 


B = 0.1244x—0.0162 


C= 1.003—0.125« 


D = 1.089% + 0.5123 


where Ep is a correction term for the earth-sun distance. 


3.1.1.7. Gueymard model. Gueymard [47,48] developed a new 
model for solar beam spectral irradiance. This model uses the 
following two possible « values: o for 2 < 0.5 um and a for other 
measurements. The f parameter is independent of the 
wavelength and is derived from the following relation: 


Be =0.5” |t5 (14) 


where Ts is the aerosol optical thickness at 0.5 um. Consequently, 
only a single measurement of the beam spectral irradiance at this 
wavelength is necessary. 


3.1.1.8. Pinazo et al. method. Pinazo et al. [40] presented a new 
method (referred to as Method II) for the calculation of the 
Angstrom turbidity coefficient based on the ratio of direct solar 
radiation to global solar radiation on a horizontal surface and on 
the “C” model of Iqbal. Method II permits the calculation of f 
without employing certain variables, such as the thicknesses of 
the precipitable water and of the ozone layer. Furthermore, 
measurements on dry temperature and relative humidity are 
not necessary. According to Method II, p can be expressed as 
follows: 


B jy P—0-12445%+0.0162)/(1.003—0.1250) 
Pin — 


m, (1.089% + 0.5123) 


The definition of D can be found in the original paper of the 
author [40]. 


(15) 


3.1.1.9. Gueymard and Vignola method. Gueymard and Vignola 
[29] proposed a semi-physical method to evaluate turbidity 
from broadband irradiance measurements and from other 
atmospheric parameters. This method should be under cloudless 
skies and requires very accurate measurements of global (or 
diffused) and direct radiations. From the number of parametric 
runs performed with SMARTS2, Kap (Kap=la/Ibn; the ratio between 
diffused and direct beam normal irradiance) can be accurately 
parameterized as follows: 


Ka = (a0 +018 +026") /(1 +436") (16) 


where the coefficients a; are functions of the zenith angle, 
pressure ozone amount, and precipitable water. 

If an experimental value of Kgp=Ig/Ipn is obtained from the 
measurements, Eq. (16) can easily be solved to obtain f as 
follows: 


p= 0.5{ [aj —4(a2—a3K ap)(do—Kab)| ° a } a2 a3Kab) (17) 


3.1.1.10. Malik method. Malik [7] modified the Pinazo et al. model 
[40] to introduce a new variable, called the “horizontal visibility,” 
to more accurately estimate the turbidity coefficient of Angstrém 
f from the ratio of horizontal direct solar irradiance to the global 
horizontal irradiance and horizontal visibility. 


r pn 2—0-12445% +0.9162)/(1.003—0.1259 
= m, (1.089% + 0.5123) 


(18) 


3.1.1.11. Utrillas et al. method. Utrillas et al. [19] developed a new 
method aimed at determining the Angstrém turbidity coefficients 
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by measuring broadband direct irradiance by using filters. 
Int,a, = Inpa Inge (19) 


where T,,, is the aerosol optical thickness at a specified spectral 
band and J/g, is the effective wavelength that corresponds to each 
filter. 


3.2. Linke’s turbidity factor (T,) 


Linke [58], [59] introduced the Linke turbidity factor (T,) and 
showed the number of dry atmospheres necessary to cause the 
attenuation of extra-terrestrial radiation produced by the real 
atmosphere [49]. T; refers to the entire spectrum or to the overall 
spectrally integrated attenuation [50]. Gueymard [51] suggested 
that the T; coefficient is not a pure turbidity coefficient because it 
also incorporates water vapor and NO2 optical depths. Never- 
theless, the T, coefficient may be useful when no data on 
precipitable water are available, making the determination of other 
coefficients difficult ([52]reported). Although T; does not explicitly 
show the effects of aerosols and water vapor, its easy measurement 
made it quite popular in meteorological stations [53]. 

T, is a key input for several models that assess the down- 
welling irradiance under clear skies. These models are used by 
several communities in the fields of renewable energies, clima- 
tology, agro-meteorology, and atmospheric pollution [54]. T; is 
also useful for the prediction of the availability of solar radiation 
and daylight illuminance under cloudless skies [55]. Although T: 
is a useful parameter in many applications, it has a serious 
drawback. T, varies with air mass even when atmospheric condi- 
tions remain constant [39]. T; normally varies from 1 to 10 [31]. In 
a pure Rayleigh atmosphere, T,=1. The value closest to this ideal 
value is achieved in extremely clear and cold air at high latitudes, 
T,=2. In a polluted atmosphere, T; can increase to 8 [56]. The 
attenuation of the radiation under a clear sky atmosphere 
becomes larger with increasing Ti. 


3.2.1. Techniques for determining the Linke’s turbidity factor (T) 

T, can be directly obtained from the pyrheliometric measure- 
ments of the beam irradiance during periods of very clear sky with 
an air mass of 2. Since this kind of experimental data is rarely 
available, T; generally becomes an estimated parameter. In addition, 
the time series of the radiation data is generally too short to allow 
estimation on a daily basis. Fortunately, long-term monthly average 
values Tım are sufficient for most applications. Several methods for 
estimating Tım values can be found in literature [27,57]. 


3.2.1.1. Linke. Linke [58] defined the integral Rayleigh optical 
thickness 6pg(m;) as the integrated optical thickness of the 
terrestrial atmosphere that is free of clouds, water vapor, and 
aerosols. He calculated 6g(m,) from theoretical assumptions and 
validated it in a very pure and dry mountain atmosphere by using 
the following formulation: 


Or(M,) = 0.128—0.054 logm, (20) 


3.2.1.2. Feussner and Dubois. Feussner and Dubois [60] published 
a series of spectral data tables that enables the calculation of 
ôr(m,), where both molecular scattering and the absorption by 
the stratospheric ozone layer are taken into account. 


3.2.1.3. Valko. Valko [61] suggested empirical fits of T; as a 
function of aerosol, water vapor content, and relative air mass. 


Valko’s formula valid for 5° < y < 65°can be expressed as follows: 
Tı = (B+0.54) [1.75 log(w/m,+0.1) +14.5]—5.4 (21) 
where B=1.07 (by Schiiepp). 


3.2.1.4. Dogniaux. Dogniaux [62] derived the following relation 
from extensive experimental campaigns, illustrating the observed 


variation in T, with solar height y (°), atmospheric water vapor 
content w (cm), and the Angstrém’s turbidity coefficient fp: 


85+/) ; : o o 
T= (rA - 0.1) +(16+0.22w)f for 5°<y<65° (22) 


3.2.1.5. Jaenicke and Kasten. Jaenicke and Kasten [63] developed 
the following linear relationship to determine the atmospheric 
turbidity factor from the threshold solar elevation that is 
measured by any given sunshine autograph with known threshold 
irradiance Ics: 


Tı = 0.154 In(Io/Ics)%¢+ 1.05 In(Io/Ics) for y 30° (23) 


where Io is the irradiance of the extra-terrestrial solar radiation and 
Ip =1375 W m7? is the mean sun-earth distance. 


3.2.1.6. Kasten. To facilitate the determination of T,(m,) from 
pyrheliometric measurements, Kasten [64] presented the following 
simple analytic formula for ôr(m,), known as Kasten’s pyrhelio- 
metric formula, based on numerical values published by Feussner 
and Dubois [59] and on international recommendations by Special 
Committee for the International Geophysical Year (CSAGI 1958): 


ôr(m,) = (9.4+0.9m,)! (24) 


3.2.1.7. Katz et al. Katz et al. [65] experimentally discovered that 
T,=2.0+19f, but they used the definition by Kasten [64] on the 
integrated optical thickness of the terrestrial atmosphere, which 
only includes Rayleigh scattering and ozone absorption. 


3.2.1.8. Kasten. T; at height z meters above sea level T,(z) was 
objectively defined by Kasten [66] as ([61] reported) follows: 


T,(Z) = 1+ 6p(Z)/dr(Z) (25) 


where or(z) is the relative optical thickness associated with 
Rayleigh scattering by the gaseous molecules in the atmosphere 
and ozone absorption, and ôp(Z) is the relative optical thickness 
associated with aerosol extinction and gaseous absorption other 
than the ozone in the stratosphere. 


3.2.1.9. Louche et al. Louche et al. [50] used more recent data on 
the spectral extra-terrestrial solar radiation and on the spectral 
absorption coefficients of the atmospheric gases, as compiled by 
Iqbal [67], to calculate the improved numerical values of the 
integral Rayleigh optical thickness dp(m;) according to the 
integration procedure outlined by Kasten [64]. Based on the 
new values of dr(m,), Louche et al. proposed a new analytic 
formula for the function dg(m,) and demonstrated that the new 
ôr(m,) values are definitely higher than the old CSAGI values used 
by Kasten, particularly, at small relative optical air masses mr. 
Louche et al. also proposed an approximate equation for 
converting old turbidity factors T,(m,) into new ones. 


Sp(m,) = (6.5567 +1.7513m,—0.1202m? +.0,0065m?—0.00013m*4) ~ 
(26) 


3.2.1.10. Abdelrahman et al. Abdelrahman et al. [53] conducted 
turbidity measurements in Dhahran, Saudi Arabia which has a 
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semi-arid climate. The inter-relationships between T; and f were 
found to be linear, similar with that in Avignon (France) and in 
Potsdam (West Germany). However, the regression coefficients 
were different for Dhahran. 


B =—0.1326 +.0.0433T, (27) 


3.2.1.11. Grenier et al. Grenier et al. [68] proposed a fourth-order 
polynomial for the function dr(m,) that is similar with Louche 
et al. [50] but with slightly different coefficients. They confined the 
range of validity to 1 < m, < 6 to achieve a better fit in this range. 


da(m;) = (5.4729 +3.0312m,—0.6329m? +.0,0910m?—0.00512m*) ' 
(28) 


Furthermore, Grenier et al. illustrated the dependence of T; 
with air mass and proposed a Linke turbidity index that was 
normalized at an air mass of 2, which they called T,4y2, and was 
therefore independent of air mass. They also investigated the 
dependence of T,4y2 on water vapor and on the Angstrom spectral 
turbidity factor and presented these dependencies in the form of 
various tables and polynomial fits. 


f =—0.120+0.071T ame 0.5 < w < 1.6 cm (29) 
B=—0.128+0.070T ay. 1.2<w<3.2cm (30) 
B=—0.138-+0.069T yy. w>2.8cm (31) 


3.2.1.12. Molineaux et al. Molineaux et al. [39] noted that the 
expressions by Louche and Grenier for dg(m,) became divergent 
for air mass measurements that were greater than 20 and 7. They 
adopted the coefficients of the original expression by Linke to 
take into account the absorption of permanent gases as follows: 


Or(m,) = 0.124—0.0656 logm, (32) 


3.2.1.13. Kasten. To compare and judge the two parameterization 
formulae for dg(m,) proposed by Louche et al. [50] and by Grenier 
et al. [68], respectively, Kasten [69] presented a new computation of 
Or(m,) according to his outlined integration procedure described in 
[64], which was repeated by Louche et al. This new computation 
disregarded several simplifications and assumptions made by 
Louche et al. and by Grenier et al., while additional parameters 
were considered and some recently improved input data were used. 


Tra(mr) = Trx(mr)[1/SpalMr)] /[1/Srx(mr)] (33) 


where T,,(m,) is the new Linke turbidity factor, T(m) is the Linke 
turbidity factor determined by the algorithm of Kasten [64], drgx(m,) 
is the integral Rayleigh optical thickness calculated using the 
formula of Kasten [64], and dgq(m;) is the integral Rayleigh optical 
thickness calculated using the formula of the new spectral 
integration, which is given by: 


Spa(tm,) = (6.6296 + 1.7513m,—0.1202m? + 0.0065m?—0.00013m*) ~" 
(34) 


3.2.1.14. Kambezidis et al. Kambezidis et al. [70] derived an 
expression for the Linke turbidity parameter in the visible 
region of the solar spectrum Ty,j;. Their expression is based on 
newly identified visible transmittances for various atmospheric 
constituents obtained by spectrally integrating wavelength- 
dependent atmospheric transmittances. 


—In(pnvis) F In(IscvisEo) 
Orvis Mr 


Tivis = 6) 5) 


3.2.1.15. Cucumo et al. Cucumo et al. [71] developed a simple and 
general calculation model for the determination of the average 
daily Linke turbidity factor in places where the daily global 
irradiation on the horizontal plane is known. Their model was 
validated using experimental data for the turbidity factor T, 
obtained in Arcavacata di Rende (Cosenza) and in Casaccia 
(Rome). The model was also experimentally validated using 
experimental T, data obtained in Arcavacata di Rende and in 
Casaccia. Moreover, polynomial correlations of the turbidity 
factor had been developed for 28 Italian localities where the 
European Solar Radiation Atlas (ESRA) provided maximum 
radiation values during clear days. Using these correlations, the 
average turbidity factor on each day of the year was calculated. 
The suggested formula for estimating T; is given by 


—InUnn) + INU scEo) 
rM; 


Ti = (36) 
where I, is the direct hourly normal energy (i.e., the incident 
hourly energy on a plane perpendicular to the Sun’s rays), dp is 
the result of the formula by Louche et al. [39] corrected by Kasten 
[69], and m, is the expression by Louche et al. [39]. 


3.2.1.16. Geiger et al. Geiger et al. [72] summarized the equations 
that describe the modified European Solar Radiation Atlas (ESRA) 
model, which was initially presented by Rigollier et al. [73] and then 
modified by Remund and Page [74] to take into account the 
variation in the integral optical depth with the altitude of the site. 


3.2.1.17. Ineichen and Perez. Ineichen and Perez [75] presented a 
new formulation of turbidity that displays a relatively stable 
behavior during a specific day. The new turbidity coefficient 
agreed with that by Kasten [64] at an air mass of 2. However, 
the contribution of absorption by atmospheric gases into the 


integral optical depth was not taken into account. The 
formulation of Tıp can be expressed as follows: 

Tup = 1 + [11.110(bEolpn)]/M; (37) 
where b is given by 

b=0.664+0.163/exp(station height, meters/8000) (38) 


The new formulation has the following advantages: (1) solar 
altitude-independent and (2) matches the original Linke turbidity 
factor at the air mass of 2. Therefore, the new formulation 
remains coherent with previous studies and is based on high 
quality data that are representative of widely differing geographic 
locations, altitudes, and climates. 


3.2.1.18. Polo et al. Polo et al. [32] presented a methodology for 
estimating daily Linke turbidity factor under clear sky conditions 
at an air mass of 2 from the global horizontal irradiance data at 
solar noon and from monthly mean values of the Linke turbidity 
factor. The method utilized the ESRA clear sky model to establish 
a functional relationship between the global irradiance at solar 
noon and the final estimation of the Linke turbidity factor. 


3.3. Schtiepp’s turbidity coefficient (B) 


Schiiepp [76] introduced the turbidity parameter B, which is 
fundamentally based similarly with £, but referred to the base 10 
instead of e and to 0.5 um wavelength instead of 1 um. 


1 In(AA) 
Ba; = log. (AA) 
m 


r 


—(TRr +Tw) (39) 


where I(^2) is the direct solar flux on the Earth’s surface, I9(A4) is 
the spectral irradiance above the atmosphere, Tę is the mean 
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Rayleigh extinction, and Tw is the selective absorption coefficient 
for atmospheric gases Oz, 03, H20, COz, and so on. 


3.3.1. Techniques for determining Schtiepp’s turbidity coefficient (B) 

Katz et al. [65] identified relationships and correlations among 
various atmospheric turbidity coefficients based on the direct 
radiation measurements in five bands of the visible region, which 
were performed in a semi-rural site. Tables 2 and 3 present the 
values of a and b coefficients in the linear regression Y=a+ bx, 
which were calculated among different turbidity parameters f, Bi, 
and Bos. 


3.4. Unsworth-Monteith’s turbidity coefficient (Ty) 


Unsworth and Monteith [77] introduced a new method for the 
estimation of the turbidity of the atmosphere. The Unsworth- 
Monteith turbidity coefficient Ty represents the number of clean 
and dry atmospheres (CDAs), which exhibit the same attenuation 
as the actual atmosphere that contains aerosols and water vapor. 
Typically, Ty varies from 0 to 1. Unsworth and Monteith used 
experimental irradiance measurements of a clean atmosphere 
and created the following formulation: 


Tyna = Tpng®XP(-Tusmr) (40) 


where I, and m, are defined earlier, Iý, represents the normal 
incidence total solar radiation under a dust-free atmosphere 
(W m~?), and exp(—Tym,) represents an aerosol spectral trans- 
mittance coefficient. With the definition of the integrated total 
optical thickness, an integrated expression of Ty that is similar to 
the expression of T; can be obtained as follows: 


Ty = mee @ (41) 


3.5. Illuminance’s turbidity factor (Tą) 


Illuminance turbidity factor Tj, was introduced by Navvab et al. 
The concept of the illuminance turbidity factor is analogous to 
that of the Linke factor with Ipn, 10, and T; substituted with Iysn, Ivo, 


Table 2 
Values of a in the linear regression Y=a+bx. 


Y/X B By B2 B3 B4 Bs Bos 
B 0 —0.005 0.008 0.002 0.007 0.003 0.001 
By 0 0.019 0.018 0.008 0.010 0.001 
B2 0 0.003 0.010 0.009 0.004 
B3 (0) — 0.006 — 0.001 0.006 
Ba 0 0.002 0.005 
Bs 0 0.006 
Bos 0 
Table 3 


Values of b in the linear regression Y=a+bx. 


Y/X B Bı B2 B3 B4 Bs Bos 
B 1 0.95 1.39 0.73 0.87 0.93 1.02 
Bı 1 1.36 0.69 0.83 0.91 1.07 
B2 1 0.51 0.61 0.68 0.72 
B3 1 1.13 1.24 1.29 
B4 1 1.10 1.19 
Bs 1 1.05 
Bos 1 


and Tj. Tj is apparently a better parameter for describing atmo- 
spheric conditions because, unlike the Linke factor, it is insensi- 
tive to water vapor content, air mass, and solar altitude angle or 
time of day. Tı can be determined from the direct measurement of 
the beam normal illuminance by Chaiwiwatworakul et al. [31] 
and Navvab et al. [78]. Illuminance turbidity can be obtained as 
follows: 


_ 1 Tysn 
Ti =— (=) (7) (42) 


The values of ô; can be calculated as follows: 
_ 0.1 
~ (1+0.0045m,) 


Illuminance turbidity has a heuristic advantage as an approx- 
imate analytical expression because it can be directly derived 
from the Angstrém formula as follows: 


Ty = 1421.68 (44) 


Oil (43) 


3.6. Volz’s turbidity factor (Ty) 


Computations of Volz turbidity are performed according to the 
usual Bouguer—Lambert-Beer Law that expresses the measured 
intensity at wavelength J as follows [79]: 


Tv = 


Me tea me] R1 (45) 


mr A 
where TÈ is the Rayleigh scattering coefficient for air molecules at 
wavelength 4, and t% is the absorption coefficient for ozone at 2. 


3.7. Kastrov’s turbidity coefficient (Tx) 


Adeyefa et al. [80] reported that Kastrov [81] obtained the 
following simple formula for the integral solar radiation on the 
surface level: 


Isc 
I = = 46 
> 4Txm*) (29) 
from which the following can be obtained: 
o1 Isc 
T= ax (F- ) (47) 


where Tc is a quantitative characteristic of atmospheric transparency. 
3.8. Rational’s turbidity factor (T,) 


Majumda et al. [82] defined a new measure of total atmo- 
spheric turbidity, termed as the rational turbidity factor, T,, to 
overcome the limitations of the Linke turbidity factor T;,. T, is 
composed of three components, namely, (i) pare and dry air (the 
basic effect), (ii) precipitable water vapor w, and (iii) aerosol 
particles in the atmosphere (dust, smoke, and haze). The formula- 
tion of T, is expressed as follows: 


Te 


1 emo oai (48) 


Mr 0.072375 


Majumda et al. [82] correlated the rational turbidity factor T, 
with the Angstrém-Schiiepp turbidity coefficient B for different 
values of precipitable water w at unit air mass to eliminate the 
effect of any possible virtual variation with air mass. They 
demonstrated that the effect of B and w were inseparable because 
the scattering of radiation by aerosol particles and the absorption 
of the same by water vapor coexist within the absorption bands of 
water vapor. Thus, the two effects cannot be isolated from each 
other. This interaction becomes more pronounced with increasing 
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particle size since the negative exponent of 4 in the Ångström- 
Schtiepp equation approaches zero for large particles. 

Majumda et al. [83] also investigated the effects of varying air 
masses m, on T; in relation to B and w and presented a formula 
that connects T, with B, w, and m, as follows: 


1—0.2347(0.869)m, 
aean je 


f= eS yey) ae aaa a! 


(49) 


4. Conclusion 
The results of this review indicate the following: 


1. In the last century, a number of atmospheric turbidity indices 
were introduced and several techniques were developed to 
determine their values. The most currently used indices 
include the Ångström turbidity coefficient £, the Linke turbid- 
ity factor T,, the Schtiepp turbidity coefficient B, and the 
Unsworth-Monteith Ty. 

2. The Angstrom exponent « is very popular because the simpli- 
city of its equation enables the interpolation or extrapolation 
of aerosol optical properties and because of its connection to 
particle microphysics. In addition, the Ångström exponent « 
also indicates the amount of aerosols and represents the 
combined effects of scattering and absorption caused by 
aerosols. 

3. The Ångström turbidity coefficient 6 has long been used as a 
turbidity index and it can be determined by different methods 
for spectral and broadband radiation measurements in large 
spectrum bands that are not affected by water vapor absorp- 
tion (280 nm < À > 630 nm). 

4. P cannot be evaluated unless « is known a priori. This 
uncertainty in « makes the determination of f less precise 
and discriminating. 

5. The Linke turbidity factor T; is important for the calculation of 
the beam solar and the diffused radiation in a specific locality. 
This parameter is also useful for the prediction of available 
solar radiation and daylight illuminance under cloudless skies. 
The Linke turbidity factor T, is also a key input to several 
models that assess down-welling irradiance under clear skies 
and that are used by several communities in the fields of 
renewable energies, climatology, agro-meteorology, and atmo- 
spheric pollution. 

6. The main problem with T, is that it is difficult to be accurately 
determined and it fails to distinguish the effect of water vapor 
from the absorption caused by aerosol. T; is also not truly 
independent of the optical air mass and it varies diurnally even 
under an unchanging atmosphere. Moreover, T; is subject to 
virtual daily variations, such as broadband aerosol optical 
depth 6, caused by the combination of different parasitic 
effects. Virtual variations in T; caused by the solar zenith angle 
y are smoother compared with those of ôa for y< 80°, but 
become steeper for larger zenith angles. A reference air mass 
value of 2 should be calculated to remove a part of the T;’s 
parasitic variation. 

7. The Angstrom and the Schiiepp turbidity coefficients, 6 and B 
respectively, have the same spectral definition. They both 
correspond to Tą, at 1000 nm (base e) and at 500 nm (base 
10), respectively, and thus, they are true turbidity coefficients 
that are unaffected by anything except by the aerosol total 
burden. 

8. The Unsworth-Monteith coefficient (or broadband aerosol 
optical depth) was found to be slightly dependent on both 
zenith angle and water vapor. Hence, caution should be 


observed when comparing data from different parts of a day or 
year. 
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